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Magnetic Fluctuations in a Charge-Ordered State of the 
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Magnetic properties in a charge-ordered state are examined for the extended Hubbard model at 
half-filling. Magnetic excitations, magnetic susceptibilities and a nuclear spin relaxation rate are 
calculated with taking account of fluctuations around the mean-field solution. The relevance of 
the present results to the observation in the 1:1 organic conductors, (TTM-TTP)l3, is discussed. 
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Organic conductors are one of the most suitable mate- 
rials for studying exotic electronic correlation effects in 
low dimensional systems. Quasi one-dimensional (1-D) 
quarter-filled conductors, Bechgaard salts (TMTSF) 2 X 
and their sulfur analog (TMTTF) 2 X, have been inten- 
sively studied and it has been found that the rich phases, 
e.g., spin-Peierls states, spin density wave (SDW) state, 
superconductivity and so on, are realized as a function 
of temperature, T, and (effective) pressure. U 

Recently, the IjJ organic material, (TTM-TT 
been synthesized© and studied experimentally. 
By extended Hiickel molecular orbital calculation, the 
transfer integral along the stacking direction is esti- 
mated as t c — 0.26e|Z and is much larger than the oth- 
ers, t a , tb <0.01eV.Q ) Therefore, the electronic struc- 
ture is highly 1-D along the stacking axis. Different 
from (TMTSF) 2 X and (TMTTF) 2 X, the compound has 
a half-filled band. Although a 1-D electron system at 
half-filling is in the Mott insulating state, the material 
shows metallic conduction for higher temperature and 
exhibits metal-insulator (MI) transition at T c = 120-160 
K.qEP Below this temperature region, the gap in the spin 
excitation has-heen observed in measurement of magnetic 
susceptibilityBtirLP and nucleac,spin relaxation rateJlP jki 
addition, X-ray measurementsBtJ> and NMR spectrurntP 
indicate the appearance of the charge order below T c . 
Thus, the MI transition of this material is followed by 
the charge order and the spin gap, and is different from 
the transition in the quasi 1-D quarter-filled organic ma- 
terials. 

Since nearest-neighbor repulsion plays a crucial role 
for the charge order, the 1-D extended Hubbard model 
(EHM) with on-site, U , and nearest-neighbor repulsion, 
V, is considered to be appropriate for studying the above 
properties. The ground state of 1-D 
has been investigated theoreticalli 
Recent investigations by numerically EZP and analytical 



calculations^ have shown the phase diagram on the 
plane of U and V as follows. In the case of large U 
and V, the SDW state is dominant for U ^ 2V, whereas 
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for U 2V the charge density wave (CDW) state corre- 
sponding to the charge order is realized where the charge- 
rich and poor sites appear alternatively. The transition 
between the two kinds of density waves is the first or- 
der. In the case of small U and V, in addition to the 
two kinds of density waves, the bond charge density wave 
(BCDW) state appears in the narrow region between the 
SDW and CDW states. The transition from BCDW to 
SDW or CDW is continuous. We note that the phase 
diagram closely resembles that by mean-field theoryfr 
where the first order transition occurs at U = 2V and 
the SDW (CDW) state is realized for U > 2V (U < 2V). 
To the best of our knowledge, however, the properties 
for finite temperature corresponding to observation in 
(TTM-TTP)I_3 and excitation spectra have not been in- 
vestigated. 

In the present paper, magnetic properties of the 
charge-ordered state in the 1-D EHM at half-filling are 
examined based on the mean-field treatment. Since the 
mean-field theory reproduces the phase diagram very 
well except the BCDW state, and the interchain cou- 
pling exists in the real materials, the treatment based on 
the mean-field solution seems to be appropriate for de- 
scribing the properties in (TTM-TTP)I_3 except the weak 
interaction region with U ~ 2V. Gaussian fluctuations 
around the mean-field solution are taken into account by 
using the path integral method. The magnetic excita- 
tion spectra, magnetic susceptibilities and a nuclear spin 
relaxation rate are investigated. 

We consider the 1-D EHM given by the following 
Hamiltonian, 7i = H.q + 7ii n t, where 
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Here, t is the transfer energy between the nearest- 
neighbor site, e k = — 2t cos ka with lattice con- 
stant a, and Qq = ir/a. The quantity, cj s (= 
l/VNJ2 k e~ lkx 'c\ s ), denotes the creation operator of 
the electron at the j-th site with spin s, n(q) = 



Ej,s 



"i C,- 



and -/V is the 

number of lattices (and electrons). We express the action 
corresponding to the above Hamiltonian by Grassmann 
algebras and decompose the interaction parts into bilin- 
ear forms by utilizing the Stratonovich-Hubbard trans- 
formation. The resulting action is written as, 

rP , 

S= / dr l y2 c *kJ T )( d r -M + £fc)cfe, s (r) 



+ 
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U 2V 
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<y 4N + N C ° S qa ^ p ( q > r M~«' T ) 



cos qa)p{q, r)n(-g,r) 
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<r{q,T)cr(-q,T) - —a(q,T)m(-q,T) }, (3) 



4jV" v "' / " v 2N~ 
where /3 = 1/T and p is the chemical potential. Here, 
we express v(0,t) = u + 5u(0,t), v(Q ,t) = vq„ + 
Sv(Qo,t) and v{q,r) = 5v{q,r) with < \q\ < Qo iy = 
p or a), and expand S up to the second order of 8v. 
Then the electron degrees of freedom are integrated out. 
The mean-field equations determining vq and vq q are 
obtained under the condition where the average of the 
terms including the first order fluctuation vanish. In 
case of U < 2V, the charge-ordered state, where <r = 
<tq = 0, becomes most stable, and the energy gap A = 
(U/2 - 2V)p Qo /N is determined by 

■ w ~ u ' 2 y:^{f{-E k )-f{E k )}, (4) 

k 



N 



with f(x) = l/(eP x + 1) and E k = ^e 2 +A 2 . The 
chemical potential is given by p = (U/2 + 2V)po/N with 
p = N. 

The effective action up to 
obtained as S c g = 



the second order is 



y—p.a cff> 



where 



quantity, A"(q, \w n ), is given as follows, 
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where 
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It should be noted that there is no coupling between 
the fluctuation with the wavenumbcr q and that with 
q — Qo although the lattice periodicity changes to 2a due 



to formation of the charge order. The reason for the 
absence of the above coupling will be discussed below. 

Charge and spin susceptibilities, Xpfei^n) and 
Xa{q,iu n ), and the local susceptibilities, x„(x i 
are given as follows, 

X P (q,iu n ) = ^{{5p(q,\u n )5p(-q,-iu n )) 

-}• (8) 
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U + 4V cos qa . 
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(Sv(q, iw n )5u(-q, -iw„)) = 
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We note that the spatially alternating component 
of the local susceptibilities, which is written as 
(-iy/N 2 J2 q (^(q,'^n)5^(-q + Qo,-iw n )>, vanishes 
due to {Sv(q, \u> n )5v{— q + Qq, —iu n )) = 0, i.e., the ab- 
sence of coupling between the fluctuation with q and that 
with q—Qo. Thus, the local susceptibilities in the charge- 
rich site are the same as those in the poor site. It is un- 
derstood in terms of the particle-hole (p-h) transforma- 
tion as follows. The Hamiltonian, cqs.(l) and (2), are in- 
variant under the p-h transformation, Cj i(r — > (— lp'ct^. 
On the other hand, the amount of charge in the charge- 
rich site changes into that in the poor site and vice versa 
by the p-h transformation. Thus, the p-h transformation 
exchanges the charge-rich site for the poor site. There- 
fore, the local quantities at the charge-rich site are equal 
to those at the poor site. In the following, we discuss the 
spin degree of freedom in detail. 

At first, we investigate the spin excitation spectra, 
which are given by the solutions of x (T (g,a)) _1 = 0. We 
show Lo(q) at T = for U/t = 3.0 and V/t = 2.0 in Fig.l 
where the solid curve represents the collective mode and 
the dotted curves, w™ m and w™ ax , express the minimum 
and maximum values of the continuum spectra, respec- 
tively. The collective mode appears just below cj™ m 
and has the maximum value for q = q m ^ir/(2a). In 
addition, the collective mode goes into the continuum 
spectra at q^O and q > q m . The locations where the 
mode touches the continuum spectra are indicated by 
the arrows in Fig.l. 

Next we discuss the magnetic susceptibility and nu- 
clear spin relaxation rate, , given by (T^T)^ 1 oc 
lim^^o I m Xa{ x ii x i'i w)/ui. As was discussed above, both 
the local spin susceptibility and (T\T)^ 1 arc indepen- 
dent of the site. These results indicate that the split of 
NMR spectra is not a Knite shift but due to the other 
origin such as a chemical shift, and the nuclear spins 
in crystallographically equivalent atoms decay with one 
kind of relaxation time. Note that it is not obvious 
whether the same conclusions are obtained or not for 
the charge-ordered states coexisting with SDW found in 
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Fig. 1. The spin excitation spectra, u>, as a function of q at 
T/t = in the case of U/t = 3.0 and V/t = 2.0. Here the solid 
curve represents the collective mode and the dotted curves, o;™ m 
and u;™ ax , express the minimum and maximum value of the con- 
tinuum spectra, respectively. The arrows indicate the locations 
where the collective mode touches the continuum spectra. 



Fig. 2. The magnetic susceptibilities, x<r(0,T) and Xv{Qo>T), 
normalized by xo(0,0) and R(T) = (TiT) - 1 normalized by 
-Rq(O) as a function of T in the case of U/t = 3.0 and V/t = 2.0. 
Here xo(0,0) and i?o(0) are respectively the susceptibility at 
q = and T = and (TiT)" 1 at T = in the absence of the 
interaction. 



quasi 1-D quarter-filling organic conductors. BEP This 
is because the p-h transformation does not exchange the 
amount of charge or spin in the rich site for that in 
the poor site though the Hamiltonian is invariant un- 
der the transformation. We show the uniform suscepti- 
bility, x<?(Q,T), the staggered susceptibility, x<x(Qo,T), 
and R{T) = (TiT)" 1 as a function of T in Fig.2. Here 
Xo(0, 0) and Ro(0) are the susceptibility and (TiT) -1 
at T = in the absence of interaction. For T > T c 
with T c being the critical temperature of the CDW, these 
quantities increase with decreasing temperature. Since 
T c is of the order of Fermi energy, 2t, Curie-like behav- 
ior in x<t(0,T) is observed. On the other hand, below 
T c , the susceptibilities decrease rapidly. The quantity, 
R(T), is enhanced just below T c and shows a rapid de- 
crease for T « T c . The behavior is the same, as Hebel- 
Schlichter peak in s-wave superconductivity^ We note 
that both quantities, x<r(Qo, T) and R(T), tend to infin- 
ity at T = T c for U = 2V, i.e., the boundary between 
the CDW and the SDW states because the transition 
between the two density waves is the first order. 

Here we consider the magnitude of the interactions. 
From the reflectance spectra, the on-site interaction, U, 
and the hopping integral, L are estimated as U ~ 0.57eV 
and t ~ 0.16eV at lOOKEP where the nearest-neighbor 
repulsion is neglected. Since there is a possibility that 
the magnitude of interaction estimated by the optical 
measurement is too big for discussing the thermody- 
namic properties f$ we may expect smaller values of 
interactions. In fact, as is shown in Fig. 3, when we 
use U/t = 0.8 and V/t — 0.6, the transition temper- 
ature is close to that for (TTM-TTP)I 3 and the mag- 
netic susceptibility for T>T C does not show Curie-like 



behavior as was observed in the experiments. In ad- 
dition, the above choice leads to the resistivity inde- 
pendent of T for hiah temperature, which is observed 
in the experimentsJjB 1 because the resistivity at half- 
filling above the charge gap behaves as T 4Kp ~ 3 f$ and 
K p = y/{\ - V/(2nt)} / {1 + {JJ + 5V)/(2irt)} ~ 3/4 in 
this case (see inset of Fig. 3). The magnitude of the in- 
teractions we have pointed out are in the weak coupling 
regime where the BCDW state appears for U ~ 2V. 
However, the estimation is considered to be sensible as 
long as the values of the interactions do not belong to the 
BCEjW region in the phase diagram determined numeri- 
callyEJ and analytically^ In fact, the values, U/t = 0.8 
and V/t = 0.6, are in the CDW region of the phase dia- 
gram. 

In conclusion, we investigated the magnetic properties 
of 1-D EHM at half-filling with taking account of fluc- 
tuation around the mean-field solution. The magnetic 
excitation spectra at T = 0, temperature dependence of 
the magnetic susceptibilities and that of the nuclear spin 
relaxation rate are discussed. The local magnetic sus- 
ceptibility and the nuclear spin relaxation rate are found 
to be independent of the site even in the charge-ordered 
state. In addition, wc pointed out the possibility that the 
magnitude of the interaction of the 1:1 material, (TTM- 
TTP)i3, is small compared to that estimated from the 
reflectance spectra. 
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